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Abstract —In substations, the presence of random transient 
impulsive interference sources makes noise highly non-Gaussian. 
In this paper, the primary interest is to provide a general model 
for wireless channel in presence of these transient impulsive 
noise for space-time signal processing problems. We assume 
a superposition of independent interference sources randomly 
distributed in space-time in a Poisson field of interferers. By 
using stochastic geometry approach, first order and second order 
statistics can be derived from basic waveforms of impulsive 
interferers. We use discrete-time series model to simulate the 
random transient impulsive waveforms. It is demonstrated that 
the amplitude distribution and density of the proposed model 
converges to a-stable distributions and their power spectral 
densities are ~ !/(/ — /o)^ where /o > 0 is a resonant frequency 
and k > 0. Measurements and computer simulations are provided 
where impulsive noise are to demonstrate the efficiency of the 
analysis. 

Index Terms —Transient impulsive noise, non-Gaussian noise 
process. Discrete-time series, Poisson field Interference, Stochastic 
geometry. 

I. Introduction 

T he presence of interferences in most environments in¬ 
valid the Gaussian noise model. For example, in the 
presence of impulsive interferences, the noise model is non- 
Gaussian ffl-El. Their impact on communication systems 
performances can be severely degraded 13, 14), IS). This 
paper is particularly focused on impulsive noise in substation 
environments. They include electrical breakdown discharges 
phenomena such as partial discharges (PD), electrical arcs 
in addition to background noise. The induced radiations are 
transient impulsive waveforms. They can occupy a wide 
frequency range which interfere with conventional wireless 
communication systems ©-ID. 

Impulsive noise modelling is an active research for the 
design of robust receivers in these environments. One of the 
most commonly used is Middleton class A model 0, a, 
©, ©, ii. It is a Poisson-Gaussian noise process where 
independent emissions of these impulses are assumed. In 121, 
0 Middleton classifies transient impulsive noise in the Class 
B model. However, it may be challenging to simulate the noise 
model due to the complex form of the amplitude distribution 
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and density which considers six parameters. An alternative 
approach can be proposed by using the well known a-stable 
distributions which is an approximation of the Middleton Class 
B 111, ITOl . In IfTTI . lfT3 . authors provide simulation methods 
to produce stable random variables. These models give suitable 
approximation of first order statistics due to strong amplitude 
and rare events of impulsive noises. In practice, the simulation 
results cannot produce transient effects induced by electrical 
discharges in substations. Thus, they are limited in terms of 
second order statistics. On the other hand, Markov-chains has 
been investigated for impulsive noise modelling in 0, ifTSl . 
States of the Markov-chain take into account various physical 
phenomena such as inter-arrival time, duration of impulses 
and their amplitude distributions. Nevertheless, the definition 
of number of states can be challenging and first and second 
order statistics cannot be written explicitly. 

In space-time signal processing problems, it is useful to 
consider the space-time distribution of interference sources 
for interference mitigation techniques and for communication 
theory in non-Gaussian noise 0, 113, US. In this paper, 
we develop a new approach to achieve a non-Gaussian noise 
model in presence of impulsive interference sources where first 
order and second order statistics can be derived and written 
explicitly. We develop an original approach for random tran¬ 
sient impulsive noise waveform modelling based on discrete¬ 
time series. Then, by using stochastic geometry approach, the 
first order of characteristic function of the random Poisson 
field of interferers can be derived. It is widely used in 
wireless communication for co-channel interference modelling 
M, ini, random access systems ESI where statistics can 
be derived. Inspired by prior works published in ED-Iisa, 
we derive first and second order statistics respectively from 
Campbell’s and Carson’s theorems in tractable forms. 

The paper is organized as follows. In section In] we provide 
a mathematical formulation of multiple interference sources in 
substation environments. The general waveform of impulsive 
noise is discussed and how the non-Gaussian noise process 
can be formulated in terms of first and second order statistics. 
In section Uni basic impulsive waveforms model is specified to 
derive the first order of characteristic function. Thus, we define 
a general impulsive waveform based on discrete-time series 
model. In section UV] the non-Gaussian process is described 
in terms of first order statistics e.g. moments-cumulants, 
amplitude distributions and densities. It is also described in 
terms of second order statistics e.g. power spectral densities. 
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Finally, in section |V] computer simulations and measurements 
are provided to demonstrate the efficiency of the analysis over 
computer simulations and measurements results. 


impulsive interference sources in substation environments. The 
first order statistics of the interfering sources can be extended 
by considering an additive background noise. 


II. Mathematical formulation of multiple 

INTERFERENCE SOURCES IN SUBSTATIONS 
A. Basic Poisson field of interferers 

Stochastic geometry approach is used to derive the first 
order of the characteristic function of interfering sources 
in presence of impulsive noise. Mathematical notations are 
inspired by Il24l . First, we need to define the statistics of 
the interference field. Let n(r, t) be a random space-time 
field as a linear superposition of individual fields randomly 
produced by activated sources (emission of radiations) in a 
domain of sources. We postulate the Poisson random field 
n(r, t) = where (5r,i is the Dirac measure on a finite 

space-time domain A C x TZ where the three dimensional 
space is considered. We note that = {r, t}^ is a set of points 
representing active sources in the space-time domain A. Thus, 
the interference field can be written as : 

n(Li) = X! ’*'utM(r.ro) (1) 

where T € £^(17, P) is an ensemble of measurable random 
waveforms of emitting sources, d(r, Tq) is a distance parameter 
where Tq is the position of the antenna and t\d{Y, Tq) is related 
to the delay of propagation induced by the position of an 
interference source and the receiver. We now assume that 
these impulsive noise are separable functions such that the 
interference field can be written as : 

n(r,t)= E E (2) 

rgAr tGWt 

where N-^ is the point process related to emitting sources and 
Nt is the point process of radiations in time domain. The 
intensity of measure of the point process Z [B) = E [\I7(i?)], 
where i? is a Borel set, has a density A('0) = A(r, t). By using 
Campbell’s theorem, the mean interference field is : 

^[n(r,f)]= [ '^r,t\d{r,ro)H'4’)d'4’ (3) 

JA 

By assuming the ergodicity of the ensemble Y such that ; 

{|d(r,ro)) ^rUt\d{r^rQ) 

where Ur.t is the basic waveform of interference sources. From 
the Laplace functional of the equation Q, the first order of 
characteristic function Q{j^) of the superposition of these 
emitting radiations is given by : 


QUO = exp y {1 - exp [-j^Ur^t]} A(V’)#^ (5) 

The first order statistics such as moments, cumulants, 
amplitude distribution and density can be derived from the 
characteristic function which depends on the definition of the 
basic waveform of interference sources. It will be specified for 


B. Interference sources in substation 

In substation environments, radiations from interferers re¬ 
ceived at the antenna are impulsive and are caused by par¬ 
tial discharges mainly in air. They can be located in HV 
equipments when physical conditions are reached to discharge 
such as presence of defects, high electric field, free electrons 
etc, II25I . In presence of multiple interference sources in the 
vicinity of the antenna, a low density of space Poisson process 
Nr can be assumed where impulsive sources are randomly 
located in the three dimensional space in far-field region. For 
the activated impulsive sources, charge particles and currents 
produced by a discharge radiate impulsive electromagnetic 
radiations. The fields E and B can be obtained from retarded 
potentials A and V respectively the magnetic potential vector 
and the scalar potential by ; 


E = -VL - — (6a) 

at 

B = V X A (fib) 

The retarded vector and scalar potentials satisfying the Lorenz 
gauge condition can be written as a wave equation : 


vV 

V^A 


1 dA 


reAfr 

(7a) 

= -do E 

reAfr 

(7b) 


where p{r,t) and J(r,f) are charge density and current 
density respectively of an activated sources Si on the emitting 
element dv^ in Vj as depicted on Fig. [T] The sum 
represents the superposition of individual source randomly 
located in the vicinity of the antenna R. The speed of 
light in the medium is represented by c, £q and po are the 
permittivity and the permeability of the vacuum respectively. 
Each activated source Si radiates electromagnetic waves in 
the medium induced by retarded potentials. A receiver R can 
received these waves on the receiving element dv/j in V/j. By 
considering successive radiations from each source in the time 
domain, the Poisson process should be extended to space-time 
process. 

Solution of equations in (|7]) is given by : 


V{r,t) = 


E L 


p{r',t- |r' - r|/c) 


47reo ; 

r'eNr 


r - r 




r'eNr ■ 


r - r 


dv (8a) 

(8b) 


The interference sources are independent such that radia¬ 
tions at the antenna surface is a superposition of independent 
impulsive noise. Generally, the receiver has a directional 
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radiation pattern not necessarily omnidirectional. The antenna 
has an effective length related to the induced voltage or current 
at the terminals to the incident field E. The receiver may also 
have RP and IF (intermediate frequency) stages e.g, low noise 
ampliher and linear hlters. As a result, the superposition of 
these radiations gives the typical waveform obtained from the 
receiver R by ; 


^t= [ aR{0X,t) *E{0X,t)dvR 
teNt 


(9) 


where aii{0,C,,t) is the aperture weighting function in sphe¬ 
rical coordinate system {r,0,Q respectively represented by 
the radial distance, the polar and the azimuthal angles. The 
aperture weighting function includes both the radiation pattern 
of the antenna and the linear impulse response of hlters. It 
can be seen as an impulse response of the receiver where the 
receiving held is converted into a time waveform alone. The 
convolution product operates for temporal impulse response. 
The antenna receives the electric held E, induced by the 
activated interferers, on the receiving element dvR in vr. 
The resulting waveform R is a superposition of independent 
impulsive noise Ut produced by activated interference sources. 
The process R is excited by a Poisson process Nt related to 
the number of impulses in time domain. It is denoted as a 
shot-noise process mi, EH. The typical impulsive noise Ut 
after any RF and IF stages of (linear) hltering is written as : 


receiver is distorted due to multipath propagation. Thus, the 
resulting impulsive noise is a random process where amplitude 
envelope and instantaneous phase are random processes. The 
propagation law may also induce the randomness of the 
amplitude scale factor. In addition, a background noiseshould 
be considered as combination of multiple independent in¬ 
terference sources below impulsive interference sources e.g, 
ambient noise from substations, thermal noise from receiver, 
etc. 

We then fully write the random process Xt as a combination 
of the shot-noise process R with an additive background noise 
rit such that : 

Xt = It Tit ( 11 ) 


C. Non-Gaussian noise process 

A common receiver design operates at a given carrier or 
center frequency. Therefore, the noise process has a resonant 
frequency such that impulsive noise is a transient signal with 
damped oscillation (see Fig. |2]). It is seen that impulsive noise 
is distorted randomly due to constructive and destructive waves 
induced by the multipath channel related to the geometry of 
interference source and the receiver. It is argued that Ut is 
generally non-stationary process where noise samples are non- 
i.i.d. 




u{0X) 




( 10 ) 


where ||m(0 , C)/r|| is the amplitude scale factor induced by 
geometry of interfering source and the receiver. u{t) is the 
amplitude envelope and (p{t) the instantaneous phase of the 
impulsive interference. In practice, receiving the signal at the 


The receiver may have a local oscillator to recover any 
desired signals in baseband. In this condition, signals can be 
demodulated at the desired resonant frequency ujq. In complex 
domain, the baseband representation is given by : 


Xt = 


(I 


Rle^iviit) 






( 12 ) 
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where the instantaneous phase of any analytic signal is ex¬ 
pressed as ip{t) = where the resonant frequency loq 

exist in (p{t). Baseband signals may be more tractable for 
impulsive noise signal processing. Indeed, the power spectrum 
density can be estimated by using classical parametric spectral 
density estimation such as Yule-Walker method. As a result, 
these impulsive noise can be reproduced by using discrete-time 
series models such as autoregressive process. 

Fig. [ 3 ] is an example of typical impulsive waveform and 
psd measured in a 735 kV substation in baseband, demod¬ 
ulated at /o = 800 MHz. It is seen that a second order of 
AR process model gives suitable estimation of the decay of 
~ 1//^ of the noise process. The non-i.i.d of noise samples 
in presence of an impulsive noise is induced by the decay 
of the power spectral density. The innovation process should 
be defined to compute distortions. The determination of first 
and second order statistics of the non-Gaussian process Xt 
strongly depends on the specification of the impulsive shapes 
Ut, see the characteristic function in the equation Q. The 
basic waveform of the impulsive noise should take account 
physical parameters such as the duration of radiations, the non¬ 
stationary behaviour of the impulsive noise Ut in which the 
amplitudes of the random process are non-i.i.d. 



(a) Impulsive noise in baseband 



Fig. 3. Impulsive noise in baseband measured in a 735 kV substation 


III. A GENERAL IMPULSIVE NOISE WAVEFORM MODEL 
USING DISCRETE-TIME SERIES 
In this section, a basic waveform of impulsive interference 
is specified. The impulse waveform at the receiver may depend 
on RF and IF stages at the receiver where the resulting 
waveform has damped oscillations generally at the carrier 
frequency. The general impulsive noise waveform model help 
to make some simple assumptions to derive first and second 
order statistics of the non-Gaussian process Xt. 

Discrete-time series models can compute these typical ran¬ 
dom waveforms observed from experimentations, i.e transient 
impulsive noise waveform with damped oscillation, damped 
exponential or a mixture of damped exponential oscillation. 
The amplitude at the present sample denoted by Ut depends 
on amplitude at the past samples denoted by Ut-i where i > 0. 
These are weighed by coefficients which give the behaviour of 
the obtained waveform Ut. The definition of these coefficients 
should be carefully defined for the stability of the process. To 
make Ut as a random process, the innovation process have to 
be a random variable to be defined. 


A. An autoregressive process for impulsive noise waveform 
modelling 

We consider real-valued random process Ut, the impulse 
shape received at the antenna produced by partial discharges 
as a discrete-time series such as an AR(p) process model is 
given by ; 


p 

Ut = 'y^,4’iUt-i + gt (13) 

i=l 

where amplitude at the past samples Ut-i are weighted by ft 
named AR(p) coefficients, et is the innovation process that 
leads to distortions of Ut. We assume a second order of the 
AR(2) process such that : 


Ut = fiUt-i4>2Ut-2 + (14) 

The AR coefficients fi and (j )2 will be defined to ensure the 
stability of the process, i.e all its roots from the characteristic 
function lie outside the unit circle. Thus, the stationarity con¬ 
ditions should be verified. The second order of the AR process 
model allows to determined roots and the autocorrelation 
function of the random process Ut easily. 


B. Definition of the innovation process 

The randomness of the process Ut is induced by the 
innovation process St. The latter should take account many 
random phenomena such as distortion of the impulsive shape 
and the non-stationary behaviour of the process linked by the 
duration of radiations received at the antenna. The definition 
of Et is based on physical assumptions : 

• When an interference source is activated, radiations re¬ 
ceived at the antenna is a superposition of the constructive 
and destructive impulsive waves caused by multipath 
effects. These can be seen as a wave distorted by an i.i.d 
random variable Et such that the equation (fl4l i is satisfied. 

• Reflectors in substations cause multiple delayed paths that 
obey to the propagation law. In this condition, the ampli¬ 
tude of the impulse received at the antenna is necessarily 
decaying with respect to time until it vanishes below 
to the background noise i.e, Et has a time-dependent 
parameter denoted by dt with the constraint that Et is 
a function that decay over time or samples. 

From these assumption, we can define Et as a white noise 
where the variance is a discrete-time function, i.e, heterosce- 
dastic white noise process Ell, m : 


St = ^tWt (15) 

where -dt is the time-dependent standard deviation of the 
white noise process Wt ^ J^{0, !)■ For those complex-valued 
random process Ut, one can consider complex-valued white 
noise process Wt. The discrete-time function -dt can be defined 
as a positive power-law or log-normal function. The latter 
takes account rise time and a fall time of impulsive noise. 
Hence we write -dt as : 
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^t = 


^0 

tatV^ 



{\0gt- fltf \ 

) 


(16) 


where at is related to the time decay of the impulse, /xj may 
refer to the time where the envelope of an impulse is maximum 
and may be related to the presence of the main path received 
at the antenna. t?o is a normalized parameter. It is convenient 
that these parameters should be set such that the rise time and 
the decay time of an impulse are much shorter than the sample 
size of the non-Gaussian process. 

A basic waveform of impulsive noise received at the an¬ 
tenna has been specified based on physical assumptions. The 
model produce impulsive noise waveforms where amplitude 
are distorted randomly by the innovation process as depicted 
on Fig. |4]where dt is a time-increment defining a sample. The 
lined curve is the real-valued impulsive noise process where 
amplitudes should decay with respect to time represented by 
the dashed curve. The process is non-stationary due to the 
time-dependent of the standard deviation of the innovation 
process St- 



Fig. 4. Impulse waveform of impulsive noise distorted by a random 
innovation process 


C. Stationarity conditions 

To ensure the stability of the process Ut, the stationarity 
conditions should be verified. By using the AR(2) process 
model, we use the Box-Jenkins modelling approach ll2^ . 
The non-stationary process is differenced until stationary 
is achieved. Hence, from the equation Gil, we write the 
difference-stationary process such that ; 


The roots depends on the value of the terms (pl + 4(^2- The 
process has a stationary solution if and only if : 


4>2 — fl < ^ 

(20a) 


(20b) 

\(j)2\ < 1 

(20c) 


These AR coefficients specify the behaviour of the wave¬ 
form of the impulsive noise Ut- They would help us to 
determine the problem statement for the determination of first 
and second order statistics of the non-Gaussian process Xt in 
section lYl 


D. Power spectral density of Ut 

By remembering that the innovation process is a heterosce- 
dastic white noise process, the power spectral density of Ut 
has a classical AR(2) psd form given by : 

where S'e(/) is the psd of e^. By using the equation dMT l in 
Appendix lAl we write the complete psd of Ut as : 

where cr^ is the variance of i?t. The variance of the white noise 
is a^ = 1. Depending on the roots of the AR process, the 
psd of Ut has different behaviour : 

> For real roots, i.e, the terms + 4(^2 > 0, if the 
characteristic equation has at least one real roots close 
to the unit circle, then Su{f) will have peak at / = 0 if 
ipi is positive. The psd will have peak at / = 0.5 if 
is negative. 

• For complex roots, i.e, the terms + Acj )2 < 0, if the 
roots are closed to the unit circle, a peak occurs near the 
resonant frequency at /o given by : 

E. Autocorrelation function of Ut 

The impulse shape has a complex form due to the ran¬ 
domness of the amplitude. It may useful to provide the 
autocorrelation function (ACF) of the process Ut- From the 
roots of the characteristic equation on GSl l, the equation GtI i 
is rewritten as : 


HL)Ut = St (17) 

where 4)(L) = 1 — fiL — (j> 2 L^ and L is the lag operator 
such that UtL'^ = Ut-i- Characteristic equation of the AR(2) 
process is given by : 

I - fiL - <I)2L^ = Q (18) 

The quadratic equation GSl) has two roots ri and ri where : 

_ (j)l± ^(1)1+ 4:(j)2 


Ut = {l-GiL)-\l-G2L)-^et (24) 

where Gi = l/ri and G 2 = l/r 2 remembering that ri 2 are 
the roots of the quadratic equation GSJ. The autocorrelation 
denoted by E[UtUt-k] = Pk of the process Ut following a 
closed form solutions : 


(Gi-G2)(1-I-GiG2) 


when ri f r 2 


1 + 


J \ 4>2 J 


(25) 


(19) 


Pk = 


when ri = r 2 
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The behaviour of the ACF pk depends on the nature of the 
roots of the quadratic function : 

• For real roots, with the constraint that |Gi| and IG 2 I < 1, 
the ACF Pk can be seen as mixture of damped exponen¬ 
tials or damped exponential oscillation that decay to zeros 
when k increases. 

• For complex roots, the ACF is a damped sinusoidal 
function where the explicit expression of pk is given by 

Ea : 

/ Sin 

where /o is the resonant frequency of the system and 
? = (1 - </' 2 )/(l + (1)2) ■ tan( 27 r/o). 

A general impulsive noise Ut has been specified by using 
AR process model. If the stationarity condition of the process 
is ensured, one can reproduce a complete random waveform 
impulsive noise Ut with damped oscillation at a desired 
resonant frequency by set the AR coefficient ifti and (j )2 such 
that the roots of the characteristic equation has complex roots. 
The impulsive noise can be represented in baseband where 
the psd Suif) should have a peak at / = 0. In this condition, 
the roots of the characteristic equation have to be real and 
particularly the first AR coefficient has to be positive, (pi > 0. 

We are now ready to derive the first and second order 
statistics of the non-Gaussian noise process Xt based on 
the equation © by using the basic waveform Ut of an 
impulsive interference from AR process model. However, it 
may be a non-trivial task since the impulsive noise shape has 
random amplitude. Some simplifying assumptions have to be 
considered. The second order statistics is the power spectrum 
of the process. It is given by the Carson’s theorem HD, ca. 

IV. First and second order statistics of the 
non-Gaussian process 

First and second order statistics are the first interest for 
implementation of threshold algorithms for signal detection 
and estimation. It may be difficult to provide exact analytical 
probability density function and power spectral densities when 
those results depend on the impulsive waveform of interference 
sources. In presence of impulsive interference sources in sub¬ 
station environments, many random impulses can be observed 
for a given time observation where inter-arrival time, energy 
of individual impulsive noise, and occurrences are randomly 
distributed. 

A. Problem statement 

We need to made some simplifying assumptions in terms 
of statistics from the basic Poisson field of interferers as well 
as in terms of basic waveforms of impulsive noise. 

1) The homogeneity of the Poisson field of interferers: 
Interference sources are generally detected in presence of 
HV equipments under voltage and the Poisson process in 
space domain can be homogeneous. However, most impulsive 
interferences are generated by AC voltages in substation envi¬ 
ronments. Thus, the interfering sources are activated whenever 
the electric field reaches the dielectric strength of the air. 


Hence, the Poisson process Nt in time domain may be inho¬ 
mogeneous and cyclostationary. However, the presence of the 
three phase voltages the superposition of the activated sources 
may become homogeneous Poisson process. In this condition, 
we assume the presence of interference sources driven by the 
three phase of voltages independently. Therefore, a constant 
density of the Poisson process A(i/>) = A is assumed. 

Furthermore, we assume a large time observation to have a 
non-negligible number of impulsive noise. As a result, we may 
write the non-Gaussian noise process Xt as a superposition of 
shot-noise process It produced by each individual interference 
source plus an additive background noise nt such that : 

We 

^ Kr.It\d{rj,ro) + (27) 

i=i 

where is a random amplitude scale factor induced by the 
geometry of individual interference source and the antenna. 
iVe is the number of the interference sources activated in the 
vicinity of the antenna. 

2) The basic waveform of impulsive noise Up First order 
statistics can be difficult to derived especially when Ut is a 
random process. Nevertheless, a suitable approximation can 
be obtained by finding an equivalent deterministic of these 
impulse shapes. In ll20l . Il26l suggest that the equivalent 
function of Ut can be determined by using the expected value 
E\Ut] denoted by 74 . From the equation ([TtI i. we write the 
equivalent impulsive noise function as : 

$(L)74 = E [£ 4 ] (28) 

The innovation process is a heteroscedastic white noise. There¬ 
fore, E [£ 4 ] = 0. Thus, the expected value of Ut is derived 
from the second order of the difference equation (|28] |. We see 
that the solution of this equation depends on the roots of the 
characteristic function. 

• For real roots, a solution of the equation (l28T l is ; 

74 = iT (e-“‘ - e-''*) (29) 

where b is related to the rise time, a is the fall time of 
the discharge 74 and K is a random amplitude factor. 

• For complex roots, a solution of the equation (l28l l is 
written as ; 

74 = AT (e““‘ - 6 “'’*) cos(27r/of -f ip) (30) 

where the resonant frequency of the system is given the 
equation (l23t and p an arbitrary phase. Here, h and a is 
not necessarily equal as find on the equation (| 2 ^ . 

It is convenient for impulsive waveforms that a and b 
are real numbers strictly positive. Moreover, we restrict the 
mathematical development of first and second order statistics 
for baseband impulsive noise. Nevertheless, one can follow 
the same approach if the damped oscillation in the equation 
(I 30 I 1 need to be considered. Therefore, first and second order 
statistics are derived based on the waveform from the equation 
(|2^ . One can recover impulsive noise with damped oscilla¬ 
tions by multiplying the impulse in baseband by a carrier wave 
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at the desired resonant frequency /q. We consider the basic 
waveform 7 ^ as a continuous-time function and we assume 
that the shot-noise process It and the background noise nt are 
independent random processes. 

B. Moments and cumulants 

The description of the shape of amplitude distributions 
and densities can be given, in some extent, by moments 
and cumulants. For example, the skewness is a measure the 
asymmetry of noise and the kurtosis, a measure of how outlier- 
prone the distribution is. We start by calculating the m}^ 
cumulant Km of the shot-noise process It- In equation ( |29] |. 
we assume K as a random variable which assume positive and 
negative values make distributions. Then, from the extension 
of the Campbell’s theorem, the cumulant Km is given by 

m ■ 


We assume a non-negative impulsive shape such that b 
higher than a. In this condition, for a finite value of m = 
{1, 2,3,4, 5, 6 }, the m-cumulants Km is a non-monotonic 
sequence if (Ff'") is negative for odd values of m, i.e Km 
is negative for odd values and Km is positive otherwise. The 
TO-cumulants can be linked with the m-moments of the non- 
gaussian model. The skewness xi the kurtosis X 2 are 
determined by the and the 4*^^ standardized moments of 
the non-Gaussian noise Xt respectively. They can be written 
in terms of cumulants such that ; 


K3 

3/2 

Kf 

(35a) 

K4 

(35b) 

^2 






= A(iT™) [ 
Jn+ 


(31) 


where Qi{s) is the first order moment generating function of 
the shot-noise process. (•) is the expectation taken over the 
distribution of the random variable. From binomial formula, 
we write the rn}^ cumulant Km as : 


• The skewness value only depends on (K^)- Indeed, if 
the skewness is |xi| < 0 , then the probability density 
function of amplitude (pdf) can be left-skewed or right- 
skewed, i.e presence of longer tail on the left or on the 
right. The pdf can be also symmetric if {K^) = 0. 
In practice, the asymmetry may be induced by random 
distortions of impulse shapes. 

• In presence of impulsive noise, if {K'^) > 0 and K 4 > 
then the kurtosis is always X 2 > 0. Hence, the pdf can 
be leptokurtic, i.e a peak around the mean and long tail 
at higher amplitude values. 


Km=X{K^) 









(- 1 )'= 


a(m — k) -\- bk 


(32) 


C. Moment generating function and characteristic function 

The moment generating function of the non-Gaussian noise 
is expressed by : 


A series can be identified with a binomial sequence. By using 
the ratio test, we prove that Km is convergent for infinite 
value of m as seen in appendix iBl Hence, Km is necessarily 
finite. We extend to the non-Gaussian noise and by assuming 
the independence between the shot-noise process and the 
background Gaussian noise. Hence, we have ; 


5s' 

grr. 


■log [Qx(s)] 


s=0 

5” 


5s' 


■ log [Q/(s)] + ^ log [(9„(s)] 


(33) 


s=0 


Finally, for an additive background noise in which an i.i.d 
Gaussian noise of zero mean and variance cr^ is assumed, 
rit ^ A/’(0, cr^). Thus, we have ; 


Qxis) = E[e-^^] 

= E 

— Qr('S)t5n('S) 


where s G TZ and Qi{s) and Qn{s) are the moment generating 
function of the shot-noise f and the additive background Gau¬ 
ssian noise nt. The moment generating function of Gaussian 
noise Qn{s) is given by ; 


( 5 „(s) = exp 




(37) 


Ki — X 


{K) jb-a) 
ab 


K2 


k=0 


(- 1 ) 




= ^E 




m 


(-1) 


a (2 — k) -\-bk 

k {Kn 


+ 


i{m — k) -\-bk 


m>2 


The generating function of the shot-noise process is more 
(34a) difficult to obtain. Nevertheless, a closed expression form of 
the moment generating function can be provided in terms of 
(34b) cumulant and by using the series expansion. We emphasize 
that the cumulant Km is finite when m goes to infinite as we 
proved in Appendix 151 Hence, we start by the series expansion 
(34c) of the cumulant generating function of the shot-noise process 
such that : 
















log[Q/(s)] = 


m—1 


= aEE 

m—1 k—0 


i-l^S 


k „m 


{K'^) 


exp 


'n 




m 


a(m — k) + bk 

(38) 


/.(.) = 

(42) 

Extended to the non-Gaussian noise, and by using the 
characteristic function of Qx{jO based on the equation ( l39b . 
we have ; 


Extended to the non-Gaussian noise, we link all values of the 
cumulant Km as defined in equations (l34l) to the moments 
of the distribution. The cumulant generating function can be 
linked by the moment generating function Qx{s) where the 
cumulant generating function is the logarithm of the moment 
generating function : 


= (43) 

The pdf fx{x) should be bounded in [0,1] Va: G 72.. We may 
rewrite the pdf by setting : 


Qx(s) 


// 

I _j_ f^mo 

^ ml 

m—1 


exp 



^ cm 

rum^ 

ml 


(39) 


where fim are moments of the distribution of the non-Gaussian 
noise such that ; 


E [Xt\ — /^i — ^1 (40a) 

E [Xf] = ^2 = K2+ MiKi (40b) 

E [xf] = P3 = K3 -I- 2K2^JLl + KifJ ,2 (40c) 

(40d) 

The cumulant generating function can be seen as an entire 
function for complex values of s where the function converges 
everywhere in the complex plane, (see Appendix I bTi. Thus, the 
characteristic function Qx{jC) can be expressed by replacing 
s = in the equation ( |39] |. 


cr^ = K2 

X — Kl 

V = - 

(7 

Erom the equation ( |43] ), we rewrite the pdf as : 


(44a) 

(44b) 





p2 


exp 



ml j 




(45) 

The expression of the pdf is complex to derived. An approxi¬ 
mation approach can help to achieve tractable forms. 

2) Series approximation of pdf: The complex form of pdf 
in equation (l43T l may be approximated by series approximation 
ifTOll . Il3n . Erom the equation (l42l) it is convenient to define ; 


Hm = (46) 

1 J-T/2 

where H{j^) is seen as the characteristic function of 7 *. In 
this condition, the pdf of the shot-noise process in the equation 
dH can be written as : 


D. Amplitude probability distribution and density of non- 
Gaussian noise process 

1) General expression of pdf: The general expression of 
probability density function of amplitude of Xt can be given 
by the convolution product of the shot-noise process It and 
the background Gaussian noise rit due to the independence of 
It and Tit- We write the pdf fxix) as : 


//(^) = [J^ exp {XTHUO - at - j^x) (47) 

Thus, from the equation (l431 l. we consider that ; 

^ = (-l)-a-i0'"(:/) (48) 

where : 


fx{x) = fl+n{x) 

= / fl{u)fn{x - u)du 
J-R 

Jn 


fxm 

0M(:y) = (27r)-l/2^g-. /2 (49) 

w V ; Q^m 

(41) Hence, the pdf of instantaneous amplitude of Xt is given 
asymptotically by collecting terms according to power of 

A- 1/2 d, El, El : 


The pdf and the characteristic function of the background 
Gaussian noise rit is well known but pdf of the shot-noise 
process It has complex form. A general form of the pdf can 
be obtained by the inverse of the Eourier transform of the 
characteristic function of the shot-noise process |[I91,E1 such 
that : 


fx{x) 


.-AT 


E 


(AT)" 


-4 




A€4Cr 


-5 


4! 


-0W(r 


2..-7 


KSCT 


72 


-0(®l(z. 


... 


(50) 
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where the first term is o(A“^/^) which is the normal distribu¬ 
tion, the second term is o(A“^) and terms within brackets is 
o(A“ 3/2) 'j’jjg approximation is based on the Edgeworth se¬ 
ries. By considering only the first term and linking the standard 
deviation a with the increment m such that am = g{o',rn) a 
function of a and m in the equation (fSOl l. one can find the 
Middleton Class A Il22l such that : 


fx{x) 



m 


(51) 


However, the Edgeworth series expansion is often inaccurate 
in the far tail of distribution ||3^ . 

3) Convergence to a-stable distribution: The non-Gaussian 
noise process can be seen as a sum of independent processes 
where the shot noise It is written as a sum of independent 
processes Ut such that : 


Ut,l+Ut,2-\ - + Ut^m , 

Xt = - - -h nt (52) 

where dm is a sequence of positive real numbers strictly 
positive. By definition, the process Xt is stable llT5l . The 
random process Ut is impulsive noise where its distribution 
is f\u\{tt) ^ where a is the characteristic exponent. 

Erom ll2Tll . Il2^ . in absolute values, one can find that the 
distribution of basic waveforms write in equations (|2^ and 
(I30ll are ^ I 7 I “ ^ where 0 < a < 2. As a result, the 
random process Xt has a a-stable distribution such that the 
characteristic function is J 35 ], m : 


QxUO = exp{jCAi- kCr (1 - J/3sign(^)?7)} (53) 

where /r is a location parameter real value, cr > 0 is a scale 
factor, /3 is the skewness parameter where — 1 <f3<l,V = 
tan(7ra/2) if a ^ 1 and p = log |^| if a = 1. 

We emphasize that the pdf of the non-Gaussian noise is a 
fat-tailed distribution with high value of kurtosis and it can 
be also asymmetric as argued. The energy and the duration 
of impulsive noise determine, to some extent, the “fatness” 
of the tail of the distribution. These parameters increase 
the probability of amplitude values higher than its standard 
deviation. The Midlleton Class A can approach the non- 
Gaussian noise however, it may be inaccurate in the tail of 
distribution. The a-stable can provide a suitable approximation 
of the amplitude distribution for those random processes which 
admit a power law decay of ^ on distribution. The 

two approximations will be compared in the section |V] The 
tail distribution is given by ; 

Fx{x) = P{X > x) = l- f fx{u)du (54) 

J —00 


E. Second order statistics : Power spectral density of Xt 

The power spectral density of the shot-noise process It can 
be given in terms of the rate A and the Eourier transform of the 
impulse response of the associated linear filter by the Carson’s 
theorem 03, Ea. The power spectral density is given by the 


autocorrelation function of the non-Gaussian noise process Xt 
is given by llT3 . 1^ , IJTI ; 


E [XtXt+r] = E [/2] + \E [ytlt+r] + al5{T) (55) 

The Carson’s theorem allow to express the power spectral 
density with these terms if the integral of the autocorrelation 
function of 74 is finite or equivalently, the integral of the psd 
is finite ESI, EQI- We assume that (A'^) < 00 . One can find 
the integral of psd of 74 is finite, S^{f )df < 00 \/f G IZ. 
In this condition, extended to the non-Gaussian process, the 
resulting expression of the psd of Xt is given by : 

Sxif) = E <5 / + , 21 + '"n (56) 

where to = 27r/. It is seen that the psd of the non-Gaussian 
noise process has a decay of ~ 1 //^. 

In this section, the first and second order statistics can be 
derived from the basic waveform of the impulsive noise. We 
proved that we can have high value of kurtosis in which 
the distribution is leptokurtic and also be asymmetric as 
discussed. Eurthermore, we proved that amplitude distributions 
and densities of the non-Gaussian process can be approximated 
by classical non-Gaussian pdf forms such as Middleton class 
A or a-stable. The power spectral density can also be derived 
where a decay of ^ 1 //^ is observed induced by waveforms 
of impulsive noise. 

V. Results and discussion 

In this section, we discuss about the validity of our theore¬ 
tical model when the electromagnetic environment has impul¬ 
sive interferers where the resulting waveforms at a receiving 
point is a succession of independent impulsive noise. 

We will start by specifying waveforms of impulsive noise 
using discrete-time series. We will define coefficients in which 
the stationarity condition is ensured for the random impulsive 
waveform process Ut, see equation (l20l l. Thus, a non-Gaussian 
noise process Xt can be fully simulated where a succession of 
random impulsive noise is excited by Poisson point process. 
Additionally, a background noise nt below the shot-noise 
process It is considered. In this condition, the first order can be 
derived where empirical amplitude distribution and density can 
be provided. We show how classical non-Gaussian noise model 
such as Midlleton class A and a-stable amplitude distribution 
and density can be appropriated vis-a-vis the simulation results 
as well as vis-a-vis real situations in substation environments. 

A. Impulsive waveforms modelling 

In section Hn] a complete random impulsive noise can be 
computed based on discrete-time series such as AR process. 
According to the equation (fl4l l. a second order of the AR pro¬ 
cess is used. We only restrict the discussion where impulsive 
noise are in baseband, i.e, a decay of ~ 1 //^ with a peak 
at / = 0. In this condition, we restrict AR coefficients such 
that their roots is real values and the stationarity condition is 
ensured, see equations (fTsT i and (l 20 l i. The innovation process 
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induced by e* is a heteroscedastic white noise process where 
the time-dependent standard deviation is given by the equation 
(US}. Parameters of the latter must be set such that the rise 
time and the decay time of a random impulsive waveform Ut 
are much shorter than the sample size of the non-Gaussian 
process. Many random impulsive noise can be simulated as 
depicted on Fig. |5] where parameters are set as follows ; 

TABLE I 

Impulsive noise shape parameters 


AR coeff. 0 

std. of et 

<^i 

(f>2 


o-t 

1.2 

-0.3 

7.0 

2.25 




500 1000 1500 2000 


0 


0 5G0 10GO 15GG 2GGG 


0 5G0 10GO 15GG 2GGG 

Sample 

(a) Amplitude 



(b) Envelope 



(c) Power spectral densities 


Fig. 5. Example of random impulsive noise Ut 


We see that the random process can generate many im¬ 
pulsive noise with random amplitudes but the power spectral 
densities have same behaviour, i.e a decay of ^ 1//^ closed 
to a Lorentzian form. The desired rise time and the decay 
time defining the duration of the impulse can be defined by 
setting parameters of the time-dependent standard deviation of 
the heteroscedastic process. One can be able to compute many 
random impulsive waveforms with many behaviours based on 
AR process as long as the stationarity condition is ensured. 


B. First and second order statistics of non-Gaussian noise 
process 

Computer simulations and measurements are provided to 
validate the analysis. 


1) Simulation setup: We are now ready to simulate a non- 
Gaussian noise process Xt in presence of non negligible 
number of impulsive noise. First of all, the electromagnetic 
environment should be specified by assuming : 

• A homogeneous random space-time Poisson field of 
interference sources where the density is an arbitrary 
constant positive value \(f)) = A < 1. 

• All radiations from interference sources emit impulsive 
noise such that parameters set on the table |I] are satisfied. 

> The energy denoted by ||C/t|P of each impulse is ran¬ 
domly distributed induced by charges and currents of 
partial discharge sources 1^ . ll^ . We choose an ex¬ 
ponential law where the energy in average, denoted by 
(II[AtIP), is above the background noise such that the 
variance ratio between the background noise the shot- 
noise process is ; 


E[I^] 


< 1 


(57) 


In this condition, we set parameters as follows : the density 
of the Poisson field interferers A(i/i) is homogeneous and 
constant and set to A = AjAr where the average interference 
sources is Ar = 5 per unit volume and the average radiation 
emissions per source is At = 5 per sample or per unit time. 
The energy is random variable exponentially distributed where 
the average value is (||C7t|p) = 10. The variance ratio between 
the background noise the shot-noise process is F = 0.1. 

2) Measurement setup: Measurement campaign is made 
in a 735 kV substation. The measurement setup includes a 
wideband antenna (0.8 to 3 GHz), RF and IF stages such as 
high pass filter, amplifier, etc. For data acquisition, we use an 
oscilloscope to capture waveforms in presence of impulsive 
noise. The sample rate is 10 Gs/s for an observation time 
at 5 ps. Details about parameters of the environment during 
the measurement campaign and the measurement setup is 
given in ioi. The obtained waveforms contain background 
noise including wireless communications from cellular or 
communication in the ISM band. We demodulate the received 
signals at the resonant frequency, /o = 800 MHz, to obtain 
waveforms in baseband. 

Results of the non-Gaussian noise process from computer 
simulation and measurement campaign in a 735 kV substation 
is provided on Fig. |6] Impulsive noise sample above back¬ 
ground noise are produced by impulsive interference sources. 


20 

15 

10 



0 1 2 3 4 5 

Sample ^lo* 

(a) Amplitude Computer simulation 


0.6 

0.4 


^ 0.2 



^ -0-2 
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-0.6 
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Time [/xs] 

(b) Amplitude S. 735 kV (40) 


Fig. 6. Example of non-Gaussian process Xt 







































11 


3) Amplitude probability distributions and densities: Am¬ 
plitude distribution and density of non-Gaussian process are 
depicted on Fig. Q and [8] which correspond respectively to 
samples from the model and from the measurement campaign 
in a 735 kV substation. It is seen that the presence of impulsive 
noise has an influence in terms of amplitude distribution 
and density. Indeed, low probability of high amplitude can 
be observed on the tail of the distribution such that we 
have fat-tailed distributions. This is due to high amplitude of 
impulsive noise especially when those amplitude distributions 
are asymptotically power law distributions. 

Classical non-Gaussian noise distributions such as Middle- 
ton Class A and a-stable distributions have these behaviours, 
i.e, leptokurtic distributions and may exhibit an asymmetry as 
seen on Fig. [8] These two distributions are compared based 
on empirical data provided by the model and measurements. 
Parameters of these distributions has been estimated from 
empirical data based on BTI - H?! . In this condition, amplitude 
distribution and density has been superposed in order to 
discuss about the quality of the fit. 



(a) Probability density 



(b) Tail distribution 


Fig. 7. Amplitude distribution and density of non-Gaussian noise process 


From the proposed model and measurements, it can be seen 
that those distributions can fit empirical data with more or 
less accuracy. The quality of the fit is determined by the 
Kullback-Liebler. It is used to measured the divergence of the 
amplitude density (pdf). Mean square error is used to compare 
tail distributions (ccdf). Results are set on the table|IIl It is seen 
that the a-stable distribution fits well the empirical data better 
than the Middleton Class A. See the KL divergence value and 
the MSE value of a-stable compared to Middleton Class A 
whatever samples from the model and measurements on the 
table El This may be explained by the approximation based 
on the Edgeworth series expansion where only the first term 
on the equation ( fSOl l is used. Therefore, a lack of accuracy is 
observed in the far tail of distribution. The a-stable distribution 
converge to the empirical data due to the definition given in the 
equation (l52l i and by arguing that impulsive noise waveforms 
Ut in absolute value can be seen as power law distributions 

f\u\{u) - 

4) Power spectral densities: Second order statistics is 
presented on Eig. |9] The psd of the non-Gaussian noise is 
estimated and smoothed by using parametric method such as 
Bug’s method m, in to observe the decay of ^ 1//^ 
induced by transient impulsive noise. It is represented by the 




Fig. 8. Amplitude distribution and density of non-Gaussian noise measured 
in 735 kV substation 


TABLE II 

Comparison of amplitude distribution and density 



Simulation 

S. 735 kV 

Model 

a-stable 

Class A 

a-stable 

Class A 

KL 

0.0037 

0.1957 

0.0111 

0.15 

MSE 

1 

o 

4.47- 10-4 

4.16 ■ 10“® 

2.37 ■ 10-4 


red curve. On the psd obtained in the 735 kV substation, 
wireless communications and harmonics can be observed at 
1.5, 2.5, 6 GHz for example. Harmonics are caused by 
interleaving artefacts and clock feedthrough from scope. 



Normalized frequency [Hz] 

(a) Computer simulation 



Fig. 9. Power spectral density of Xt 


VI. Conclusion 

In this paper we develop a non-Gaussian noise model in 
presence of transient impulsive noise in substation environ¬ 
ments. We use Poisson field of interferers in which impulsive 
transient interference sources are space-time Poisson process. 
Based on stochastic geometry, first and second order statistics 
can be derived. In presence of impulsive noise, it is proved that 
the amplitude distribution and density can be approximated 
by classical non-Gaussian noise such as Middleton Class A 
and a-stable distributions. It is seen that the latter is a better 
approximation than the Middleton Class A due to its approxi¬ 
mation by using Edgeworth series expansion. Basic impulsive 
waveform is specified by using discrete-time series where 
the innovation process is heteroscedastic to ensure both the 
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the randomness and the transient behaviour of the impulsive 
interference sources. It is proved that the non-Gaussian noise 
can be expressed as a second order statistics where the power 
spectral density is a decay of ~ I//*- Computer simulation 
and experimental data are provided to show the validity of the 
analysis. 

In Future works we will discuss about the reproducibility of 
the model vis-a-vis the experimentations. Physical parameters 
such as number of activated interfering sources in the envi¬ 
ronment, the variance ratio between background noise and the 
shot-noise process and duration of impulsive noise have to be 
estimated from the desired environment and validated in terms 
of first and second order statistics. 


L = lim 

k^oo 


Zk+1 

Zk 


< 1 


For the binomial coefficients we have ; 


(63) 



(64a) 

(64b) 


In this condition, based on the equation (l32]) . we write L such 
that : 


Appendix A 

Power spectrum density oe heteroscedastic 

PROCESS et 

From the innovation process St defined in equation ( fTSl l, we 
start by calculating the autocorrelation function as : 


E [etet-k] = E [^t^t-kWtWt-k] (58) 


L = lim 

k—¥oo 


(—1 )(to — k) 

am + {b — a){k + 1) 

k + 1 

am + {b — a)k 


(65) 


We define b higher than a. Moreover, since k < m and positive 
integers, m is necessarily infinite, we determine that L = 0 < 
1. It is proved that Zk is convergent and thus, the series Km is 
convergent. 


The white noise process W* is Ltd such that : 

E [WtWt-k] = 0 (59) 

for all values of k ^ 0. By assuming that dt and Wt are 
independent variables, we write the equation (fSST l as ; 

E [etst-k] = E E [WtWt] 

= E a^5(0) 


B. The radius of convergence of the cumulant generating 
function 

From the power series expansion of the cumulant generating 
function is given by the equation dBSl) . where s £ C is complex, 
the radius of convergence of the power series can be discussed. 
A power series converge for some values of the variable s and 
may diverge for others. Thus, the radius of the convergence 
can be calculated from Cauchy-Hadamard theorem’s : 


where is the variance of the white noise. The Fourier 
transform of the autocorrelation function of St allow us to 
write as discrete convolution product between the variance of 
the white noise and the psd of dt denoted by Ss{f) ■ 


r 


-1 


lim sup I Km I "* 

m—^oo 


( 66 ) 


The radius of the convergence can be calculated from the 
ratio test of Kmltrd. : 


*5'e(/) — S^{f) * cr^ 

= f crwMf)df (61) 

Jn 

2 2 

where (t| is the variance of dt- 

Appendix B 

About the convergence oe the cumulant Km and 
THE cumulant GENERATING EUNCTION 

A. Convergence of cumulant Km 

We postulate that k < m where k and m are positive integer 
and m > 0. From the equation (|32] |. we denote (z^) a sequence 
of non-zero real values sequence as ; 


r ^ = lim 


= lim 

m—foo 


^m+1 


(m -I- 1 )n„ 

m+1 

E (™+i) [a{m + l-b) + bk]\ 




(^ + 1 ) E (T) [a(™ -b) + bk] I (K^) 

fc=0 

lim = 0 

^->oo 771-1-1 


(67) 


where r is the radius of convergence which is infinite, r ^ oo, 
i.e the cumulant generating function converges everywhere in 
the complex plane. Therefore, it is an entire function. 


Zk 



(- 1 )'= 


(AT’") 

a{m — k) + bk 


(62) 


By using the ratio test, the convergence of the series Km is 
ensured if and only if : 


Appendix C 

Power spectral density oe jt 

The power spectral density of 77 is given by the Wiener- 
KJiinchine theorem : 
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S'i{f)= / E[jtjt+r]e^‘^*dt 
Jn 

SM) = {K^) 

= {K^) 

{K^){b-af 
(a^ + a;^)(6^ + w^) 


1 1 
a + juj 5 + joj 
b — a 

(a + juj)(b + juj) 


( 68 ) 


where w = 27r/. The psd of jt is finite such that S^{f) < oo 
for all values of / G 7^ by assuming that < oo. As a 

result, the integral of the psd S^{f) is finite. 
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